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0.72  degree  Global  Domain 


Pan-Am  Global  Grid 

o  0.72  degree  equatorial  Mercator  78S-47N 
o  Arctic  bi-polar  patch  above  47N 

*  Low  resolution  global  had  patch  at  59N 

*  Can’t  include  Hudson  Bay 

o  Double  latitudinal  resolution  near  the  equator 
o  Halve  latitudinal  resolution  in  Antarctic 

Coastline  at  50m  isobath 

o  Closed  Bering  Strait 

o  No  Sigma  (terrain-following)  vertical  coordinate 

Same  26-layers  as  0.08  degree  Atlantic 

Two  configurations  for  fixed  coordinates: 

o  First  Z-level  3  m,  increases  1 ,125x  up  to  12  m 
o  First  Z-level  3  m,  increases  1 ,200x  up  to  60  m 

Also  a  40-layer  Z-coordinate  case: 

o  First  Z-level  1  m,  increases  1 ,190x  up  to  878  m 


47N:  SCPY  (km) 
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47N:  Grid  Aspect  Ratio  (SCPX/SCPY) 
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Three  global  resolutions 


0.72,  0.24  and  0.08  degrees 
0.72  degrees  is  inexpensive  and  non-eddying 
o  Only  resolution  run  so  far 
0.08  degrees  is  expensive  and  eddy  resolving 

o  Target  resolution  for  operational  use  at  NAVO 

o  Concentrate  on  Atlantic  and  Pacific  at  0.08 
FY04 

0.24  degrees  is  (marginally)  eddy  permitting 

o  Primary  resolution  in  FY04 

Starting  point  for  0.24  and  0.08  is  a  0.04 
coastline/bathymetry 

o  Interpolated  from  NRL  2  minute  bathymetry 

o  Extensive  quality  control  in  straits  and  near 
coastlines 


o  Still  under  preparation 


0.72  degree  Global  Standard  Configuration 

KPP  mixed  layer 

Energy-Loan  ice  model 

Sigma-theta  (some  sigma2  runs) 

Horizontal  diffusion  chosen  to  suppress  eddies 

Initialize  from  GDEM3 

ECMWF  Reanalysis  monthly  mean  forcing 

o  Plus  6-hrly  wind  anomalies  from  sep94-sep95 

Longwave  correction  w.r.t.  ECMWF  SST 

Inexpensive  approximation  to  COARE  2.6  bulk  heat 
flux  parameterization 

Monthly  means  of  15  largest  rivers  via  precip  bogus 

Strong  relaxation  to  monthly  GDEM3  SSS 

o  “30  days  in  30  m”  e-folding  time 
o  Necessary  to  prevent  SSS  drift 
o  In  addition  to  E-P  forcing  (monthly  P) 


Longwave  Radiation  and  SST 


Longwave  Radiation  is  sum  of: 

o  Upward  blackbody  longwave  radiation 

*  Qbb  =  -0.98  (5.67  x  10“8)  (Ts  +  273.16)4 

o  Downward  atmospheric  longwave  flux 

*  Highly  dependent  on  cloudiness 

*  Unknown  dependence  on  SST 
(assume  independent) 

If  longwave  was  calculated  using  a  SST  of  Tso\ 

°  Qiw(Ts)  =  Q,JTso)  +  Qbb(Ts)  —  Qbb(TS0 ) 

°  Q,  JTS)  =  Qiw(Tso)  +  Qbb  —  Tso) 

O  Q'bb  =  -0.98  (5.67  x  10“8)  4  (Ts  +  273.16)3 

Ocean  Model  Intercomparison  Project  includes 
(Qbb(Ts)  -  Qbb(Tso ))  as  a  longwave  correction 

HYCOM  uses  the  approximation  (in  W/m2)\ 

O  Q'bb  =  -4.506  -  0.0554  Ts 

This  is  similar  to  “30  days  in  3.5  m”  SST  relaxation 
o  lOx  weaker  than  typically  SSS  relaxation 


SST  Metrics  I 


Run  for  5  years  and  form  monthly  means 

o  Takes  two  days  on  64  IBM  POWER4  opus 
o  25  year  run  gives  “similar”  SST 

Compare  monthly  SST  to  Reynolds  and  Smith 
climatology 

o  Monthly  anomalies 
o  Annual  mean  difference 
o  RMS  difference 
o  Correlation  Coefficient 
o  Skill  Score 

*  Correlation  squared  -  Unconditional 
Bias  -  Conditional  Bias 

*  Maximum  is  1 ,  but  minium  is  -infinity 

*  Measure  of  error  w.r.t.  seasonal  cycle 
(i.e.  w.r.t.  standard  deviation) 

*  Use  a  minimum  of  1  degC  for  standard 
deviation 

•  Still  get  poor  skill  scores  near  equator 


SST  Metrics  II 


Purpose  of  comparison  is  to  find  “good  enough” 
configuration 

o  Assume  that  “skill”  on  climatological  forcing  is 
maintained  on  interannual  forcing 

o  Is  monthly  thermal  climatological  forcing  enough? 

o  NLOM  experience  suggests  that  this  is  OK,  but 
can’t  be  certain  until  we  run  more  interannual 
cases  with  HYCOM 

Targets: 

o  Annual  mean  error  <  0.5  degC 
o  Correlation  Coefficient  >  0.6 
o  Skill  Score  >  0.3 

Use  zonal  averages  to  reduce  amount  of  data 

o  Average  not  necessarily  best  statistic 

*  A  few  large  negative  skill  scores  can 
dominate  the  average 

o  Same  targets  as  for  full  field 


Simulation  History 


Expt  3.3: 

o  Best  as  of  February  2003 
o  Levitus  climatology 
o  Annual  rivers 

o  “Longwave”  via  SST  relaxation 

Expt  9.0: 

o  40  Z-level  case 
o  PLM  vertical  remapping 

Expt  1 1 .0: 

o  Standard  sigma-theta  case 
o  3m-12m  Z  levels 
o  Thin  deep  isopycnal  layers 


ZONAL  SST  MEAN  ERROR,  HYCOM  VS  Reynolds&Smith 
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ZONAL  SST  Correlation,  HYCOM  VS  Reynolds&Smith 
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ZONAL  SST  Skill  Score,  HYCOM  VS  Reynolds&Smith 
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03.3  vs  R&S  SST:  Mean  Error 


11.0  vs  R&S  SST:  Mean  Error 
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09.0  vs  R&S  SST:  Mean  Error 
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1 1 .0  vs  R&S  SST :  Mean  Error  yr  25 
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03.3  vs  R&S  SST:  Skill  Score 


11.0  vs  R&S  SST:  Skill  Score 
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09.0  vs  R&S  SST:  Skill  Score 
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1 1 .0  vs  R&S  SST :  Skill  Score  yr  25 


Interannual  SST  Comparisons 


Identify  year-long  time  series  at  fixed  locations 
o  Always  SST 

o  Sometimes  atmospheric  fields 
o  Sometimes  subsurface  T  and/or  S 

Compare  observations  to  free-running  and 
assimilative  ocean  models 

o  HYCOM  0.72  global,  free  running  (ECMWF 
or  NOGAPS) 

o  NLOM  1/8  near-global,  free  running  ECMWF 

o  NCOM  1/8  global,  free  running  FNMOC  and 
assimilative 

NLOM/NCOM  1/8  degree  is  1/6  degree  equatorial 
o  4x  finer  than  HYCOM 


Number  of  buoys  Number  of  buoys 


Free-running  HYCOM,  NCOM  and  NLOM  with  no  assimilation  of  SST 
Comparisons  with  86  year-long  daily  buoy  time  series  (1998-2000) 
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|  NLOM  (Median  R=0.81) 
|  NCOM  (Median  R=0.88) 
HYCOM  (Median  R=0.85) 


(-O.l^O.O)  (0. 1-0.2)  (0.3-0. 4)  (0.5-0.6)  (0.7-0.8)  (0.9-1.0) 
Correlation  coefficient 


|  NLOM  (Median  RMSD=0.95°C) 
|  NCOM  (Median  RMSD=1.54°C) 
HYCOM  (Median  RMSD=0.91°C) 


(0. 4-0.6)  (0. 8-1.0)  (1.2-1. 4)  (1.6— 1.8)  (2.0-2.2)  (2.4-2.6) 
Root  mean  square  difference  (°C) 


1/8°  NLOM  uses  wind  and  thermal  forcing  from  ECMWF 
1/8°  NCOM  uses  wind  and  thermal  forcing  from  NOGAPS 
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Error  statistics  with  respect  to  buoy  (365  days) 

ps:  HYCOM  and  NCOM  are  free — running  simulations 
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HYCOM  (KPP)  vs  Buoy  at  38W,15N  in  1998 
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Conclusions  (SST) 


Global  0.72  degree  SST  is  similar  to  that  in  higher 
resolution  Atlantic  and  Pacific  simulations  (not  shown) 

KPP  is  performing  well  in  HYCOM 

Thinner  deep  isopycnal  layers  are  a  major 
improvement  on  the  equator 

Skill  in  southern  mid-high  latitudes  and  northern 
(Atlantic)  high  latitudes  needs  improving 

Most  of  the  SST  error  is  in  the  annual  mean 

Not  yet  clear  how  much  is  due  to  forcing 
and  how  much  due  to  KPP 

o  If  it  is  due  to  forcing,  we  can  apply  a  correction 


Long  Term  Trends 


Most  cases  run  only  5  years, 
but  1 1 .0  has  run  25  years 

No  significant  trend  in  SST  or  SSS 

No  significant  trend  in  heat  flux  or  E-P 

Probably  too  little  sea  ice 

o  Antarctic  extent  shrinks  over  time 

Some  basin-wide  averages  show  a  trend 

o  SSH  increasing  by  5  mm/year, 
rate  changing  -0.1  mm/year 

*  Steric  change,  i.e.  lighter  average  density 
o  Average  T  warming  by  0.6  degC/century 

*  Even  though  net  heat  flux  is  cooling 

o  Average  S  freshening  by  0.05  psu/century 

*  Net  E-P  is  also  freshening 

We  need  to  find  a  way  to  equilibrate  SSH  (density) 


Global  0.72deg,  Basin-Wide  Mean  Heat  Flux  Global  0.72deg,  Basin-Wide  Mean  Surface  P-E  Flux 
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Global  0.72deg,  Percentage  Global  Ice  Cover  Global  0.72deg,  Basin-Wide  Mean  Temperature 
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Thermohaline  Circulation 


North  Atlantic  Overturning  StreamFunction 

o  Initially  very  strong  in  1 1 .0  (sigma-theta) 
o  Weakens  over  time 

Sigma2*  would  presumably  be  better 

o  But  reference  state  isn’t  simultaneously 
stable  in  Antarctic  and  Labrador  Sea 
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Depth  of  Maximum  of  Atlantic  Overturning  Streamfunction 
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Thermobaricity  in  Global  Models 

Need  sigma2  to  represent  AABW 

Need  thermobaricity  to  get  a  good  SSH  with  sigma2 

MICOM/HYCOM  include  thermobaricity  via  “virtual 
potential  density”  (sigma2*) 

o  Compressibility  coefficient  from  a  reference  T&S 

Scheme  is  unstable  if  actual  T&S  are  very  different 
from  reference  T&S 

No  single  reference  T&S  works  globally 
Tried  using  combination  of  two  reference  states 

o  3  degC  and  35  psu  north  of  30N 
o  0  degC  and  34  psu  south  of  30S 
o  weighted  sum  between  30N  and  30S 

Tested  the  method  by  shrinking  the  transition  zone 

o  Expt  16.1 :  30S-30N  (as  above) 
o  Expt  16.3:  30S-25S 
o  Expt  16.4:  25N-30N 

This  approach  does  not  appear  viable 

What  else  can  we  try? 


ZONAL  GLOBAL  FSD  ANOMALY,  HYCOM  16.1  vs  HYCOM 
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Future  Work 


In  FY04,  most  of  global  effort  will  be  at  0.24  degrees 
o  Perhaps  a  0.08  degree  “demo”  run 
Continue  to  concentrate  on  “free  running”  SST 
o  KPP  vs  GISS 

o  Corrections  for  atmospheric  biases 
o  Interannual  comparisons  to  buoys 

Other  issues: 

o  Add  CICE  sea  ice  model 
*  How  to  evaluate  sea  ice  skill 
o  SSH/density  drift 
o  Thermobaricity 
o  SST  assimilation 

o  What  is  the  optimal  surface  salinity  forcing 


